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25 ABSTRACT 

26 Coronavirus (CoV) non-structural protein 14 (nspl4) is a 60 kDa protein encoded by 

27 the replicase gene that is part of the replication-transcription complex. It is a 

28 bifunctional enzyme bearing 3'-5' exoribonuclease (ExoN) and guanine-N7- 

29 methyltransferase (N7-MTase) activities. ExoN hydrolyzes single- and double-stranded 

30 RNAs and is part of a proofreading system responsible for the high fidelity of CoV 

31 replication. Nspl4 N7-MTase activity is required for viral mRNA cap synthesis and 

32 prevents the recognition of viral mRNAs as “non-self’ by the host cell. In this work, a 

33 set of point mutants affecting different motifs within the ExoN domain of nspl4 was 

34 generated, using transmissible gastroenteritis virus as a model of Alphacoronavirus. 

35 Mutants lacking ExoN activity were non-viable despite being competent in both viral 

36 RNA and protein synthesis. A specific mutation within zinc finger 1 (ZF-C) led to a 

37 viable virus with growth and viral RNA synthesis kinetics similar to that of the parental 

38 virus. Mutant rTGEV-ZF-C caused decreased cytopathic effect and apoptosis compared 

39 with the wild-type virus and reduced levels of dsRNA accumulation at late times post- 

40 infection. Consequently, the mutant triggered a reduced antiviral response, which was 

41 confirmed by evaluating different stages of the dsRNA-induced antiviral pathway. The 

42 expression of IFN-J3, TNF, and interferon-stimulated genes in cells infected with mutant 

43 rTGEV-ZF-C was reduced, when compared to the parental virus. Overall, our data 

44 revealed a potential role for CoV nspl4 in modulation of the innate immune response. 

45 

46 IMPORTANCE 

47 The innate immune response is the first line of antiviral defense that culminates in the 

48 synthesis of interferon and proinflammatory cytokines to control viral replication. CoVs 

49 have evolved several mechanisms to counteract the innate immune response at different 
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50 levels, but to date the role of CoV-encoded ribonucleases in preventing activation of the 

51 dsRNA-induced antiviral response has not been described. The introduction of a 

52 mutation in zinc finger 1 of the ExoN domain of nspl4 led to a virus that induced a 

53 weak antiviral response, most likely due to the accumulation of lower levels of dsRNA 

54 in the late phases of infection. These observations allowed us to propose a novel role for 

55 CoV nspl4 ExoN activity in counteracting the antiviral response, which could serve as 

56 a novel target for the design of antiviral strategies. 
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57 INTRODUCTION 

58 Coronaviruses (CoVs) are enveloped, single-stranded, positive-sense RNA viruses 

59 belonging to the Coronaviridae family within the Nidovirales order (1). CoVs are 

60 frequently associated with respiratory and enteric diseases in humans, livestock, and 

61 companion animals (2, 3). CoVs have been divided into four genera: Alphacoronavirus, 

62 such as human coronavirus 229E (HCoV-229E) or the porcine transmissible 

63 gastroenteritis virus (TGEV); Betacoronavirus, that includes the severe acute 

64 respiratory syndrome coronavirus (SARS-CoV) and the recently emerged Middle East 

65 respiratory syndrome CoV (MERS-CoV), both causing pneumonia and having a high 

66 mortality rate; Gammacoronavirus, including infectious bronchitis virus (IBV); and 

67 Deltacoronavirus (4). CoVs contain the largest known genome among RNA viruses, 

68 consisting of a single-stranded, positive-sense, 5'-capped and polyadenylated RNA 

69 molecule of 27-31 kb in length (5). The first two-thirds of the genome contains the 

70 replicase gene, which is comprised of two overlapping open reading frames (ORFs), la 

71 and lb. After infection, viral genomic RNA (gRNA) is directly translated to yield two 

72 polyprotein precursors, ppla and pplab, the latter requiring ribosomal frameshifting 

73 near the 3' end of ORFla for its translation (6). Subsequently, the two polyproteins are 

74 cleaved by ORFla-encoded proteases to release a total of 16 nonstructural proteins 

75 (nspl to nspl6) (7, 8). These nsps assemble together with cellular factors to form a large 

76 replication-transcription complex (RTC) associated with membrane structures derived 

77 from the endoplasmic reticulum (9-11). The RTC is required for CoV RNA synthesis, 

78 which is comprised of genome replication and transcription. CoV transcription is a 

79 discontinuous process that yields a nested set of subgenomic mRNAs (sg mRNAs) that 

80 serve as templates for translation of the viral structural and accessory proteins (8, 12, 

81 13). The CoV RTC is extremely complex, and apart from the RNA-dependent RNA 
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82 polymerase (RdRp) and helicase activities common to many RNA viruses, CoVs 

83 encode a unique set of RNA-modifying activities such as the 3'-5' exoribonuclease 

84 (ExoN) and the uridylate-specific endoribonuclease (NendoU), a recently discovered 

85 nucleotidyltransferase (14), or enzymatic activities related to the synthesis of the cap 

86 structure, such as RNA 5'-triphosphatase (RTPase), N7-methyltransferase (N7-MTase), 

87 and 2'-0 methyltransferase (20-MTase) (7, 15, 16). 

88 CoV nsp!4 is a bifunctional enzyme that harbors both ExoN and N7-MTase activities 

89 (17-19). The amino-terminal part of nspl4 includes the ExoN active core, which is 

90 divided into three motifs: I (DE), II (E) and III (D). Due to this characteristic, CoV 

91 nspl4 is included in the DEDD exonuclease superfamily, which comprises cellular 

92 enzymes that catalyze DNA proofreading (20). Nspl4 ExoN has been proposed to have 

93 a critical role in CoV replication and transcription, as mutants lacking this activity 

94 showed an important reduction in viral RNA synthesis or accumulation (17, 21). Nspl4 

95 hydrolyzes single- and double- stranded RNA (ssRNA and dsRNA) to final products of 

96 8-12 nt and 5-7 nt, respectively, and its nucleolytic activity is enhanced up to 35 fold by 

97 interaction with nsplO (17, 22). In addition, nspl4 ExoN activity was proposed to be 

98 part of the RNA proofreading machinery during CoV replication (20, 23), which would 

99 be required for the replication and maintenance of the large CoV genome. In fact, only 

100 larger-sized Nidovirus genomes encode ExoN activity (24). Mutations in the ExoN 

101 active core abolishing the exonuclease activity led to CoVs that exhibited a 15-20 fold 

102 increase in replication errors in Betacoronavirus such as MHV and SARS-CoV (21, 23, 

103 25). As a proofreading component, ExoN should be involved in the removal of 

104 misincorporated nucleotides. Indeed, nspl4 activity efficiently removed mismatched 3'- 

105 end nucleotides mimicking RdRp misincorporation products (22). Moreover, mutants 

106 lacking ExoN activity showed greater sensitivity to the mutagen 5-fluorouracil, in 
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107 contrast to CoVs with ExoN activity, which are considered resistant to lethal 

108 mutagenesis (26). 

109 The carboxy-terminal part of nspl4 contains N7-MTase activity, involved in the 

110 addition of a methyl group to the cap guanosine at the N7 position, leading to formation 

111 of the cap-0 structure (18, 27). In general, this cap-0 is critical for efficient export, 

112 translation, and stability of mRNAs. In addition, the methylation of N7 would be 

113 required for the subsequent methylation at the 0-2' position, which is essential for 

114 prevention of viral RNA recognition by the host immune system (28-30). 

115 The bifunctional nspl4 is part of the large multi-subunit polymerase complex described 

116 as the core of the CoV RTC, which integrates RNA polymerization, proofreading, and 

117 cap-modifying activities into a multifunctional protein assembly (16, 31). 

118 The recent crystallization of the SARS-CoV nspl4 identified the presence of three zinc 

119 fingers (ZFs) within the nspl4 structure (32). A ZF is a small, independently folded 

120 domain that is structured around a zinc ion, which is coordinated through cysteine and 

121 histidine residues. ZFs are structurally diverse, with more than 40 types of annotated 

122 ZFs in UniProtKB, and are found in proteins that perform a broad range of functions, 

123 including essential cellular processes such as replication, signaling, cell proliferation, or 

124 apoptosis (33). ZFs usually function as interaction modules and bind to a wide variety 

125 of compounds, such as nucleic acids, proteins, and small molecules (34). While 

126 interaction of different ZF types with DNA has been well characterized in terms of 

127 affinity and specificity, ZF interactions with RNA remains poorly characterized (35). 

128 Nevertheless it has been shown that modification of a ZF sequence may alter its binding 

129 to an RNA substrate (36). 

130 Viral dsRNA produced as a replication intermediate is a pathogen-associated molecular 

131 pattern (PAMP), which mediates the activation of well characterized antiviral 
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132 mechanisms leading to shutdown of protein synthesis, stimulation of host innate 

133 immunity for initial detection of pathogens, and subsequent activation of adaptive 

134 immunity (37, 38). CoVs have evolved mechanisms impairing the activation of the 

135 innate immune response at different stages. These mechanisms include the production 

136 of several viral proteins that interfere with type I interferon (IFN) and proinflammatory 

137 cytokine production (39) and the induction of double-membrane vesicles, where dsRNA 

138 is shielded from innate immune sensing by cytoplasmic receptors (10). 

139 As a key enzyme in cap formation, CoV nspl4 was proposed as an IFN antagonist (39). 

140 In fact, overexpression studies have revealed that nspl4 acts as an IFN antagonist (40), 

141 although it has not been determined whether this activity is linked to its ExoN or N7- 

142 MTase domains. The second CoV-encoded ribonuclease, NendoU, is also a robust IFN 

143 antagonist when overexpressed (41). The possibility that CoV ribonucleases could 

144 degrade dsRNA PAMPs has been suggested as an additional mechanism preventing 

145 IFN induction during CoV infection (42). Nevertheless, the role of these proteins in 

146 preventing the innate immune response has not been studied so far. 

147 Current knowledge about CoV nspl4 is mainly derived from study of the 

148 Betacoronaviruses. In contrast, there is limited information on the role of 

149 Alphacoronavirus nspl4 during infection, as HCoV-229E mutations which abrogated 

150 ExoN activity were lethal (17). In this work we used another Alphacoronavirus, TGEV, 

151 to analyze the role of the nspl4 ExoN domain in terms of virus viability and induction 

152 of antiviral response. Our results show that mutations abolishing ExoN activity resulted 

153 in only a modest decrease in viral RNA synthesis, but recovery of infectious progeny 

154 failed, something that may be a general characteristic for alphacoronaviruses. A specific 

155 mutation in the ZF1 produced a viable virus causing a decrease in accumulation of 

156 dsRNA intermediates at late times post-infection. Consequently, this mutant virus 
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157 triggered a reduced antiviral response and apoptosis in comparison to the parental virus, 

158 indicating that CoV nspl4 plays a role in modulation of the innate immune response. 

159 

1 60 MATERIALS AND METHODS 

161 Cells. Baby hamster kidney cells (BHK-21) (ATCC CCL-10), or BHK-21 cells stably 

162 transformed with the gene coding for porcine aminopeptidase N (BHK-pAPN) (43) 

163 were grown in Dulbecco's modified Eagle's medium (DMEM) supplemented with 5% 

164 fetal calf serum (FCS). Both cell lines were transfected with the Sindbis virus replicon 

165 pSINrepl (44) expressing TGEV nucleocapsid (N) protein, obtaining BHK-N or BHK- 

166 pAPN-N cells, respectively. G418 (1.5 mg/ml) and puromycin (5 pg/ml) were used as 

167 selection agents for pAPN and pSINrepl, respectively. Recombinant TGEV viruses 

168 obtained in this work were grown in swine testis (ST) cells (45) using DMEM 

169 supplemented with 10% FCS. 

170 

171 Plasmid constructs. cDNAs of TGEV-derived replicons and infectious viruses (46, 47) 

172 were generated by PCR-directed mutagenesis. To generate nspl4-ExoI, nspl4-ExoIII, 

173 nspl4-ZF-H, nspl4-ZF-C and nspl4-N7MTase mutant sequences, two overlapping 

174 PCR fragments were obtained by using as a template the plasmid pBAC-TGEV (47), 

175 containing the full-length TGEV genome (GeneBank accesion number AJ271965). The 

176 5' fragment was obtained by using the forward primer nspl3-finVS (5- 

177 CATGTGTGATAGAACTATGTATGAGAATCTTG-3') and the specific reverse 

178 primer shown in Table 1. The 3' fragment was obtained by using the specific forward 

179 primer shown in Table 1 and, in all cases, the reverse primer nspl5prin-RS (5- 

180 CCATTATTTTGTCAGCAATAACAGCAG-3'). Full-length amplicons were generated 

181 by overlap extension of 3' and 5' PCR products using primers nspl3-finVS and 
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182 nspl5prin-RS. In all cases, full-length amplicons were digested with BstBI and Xbal 

183 restriction enzymes and cloned into the same sites of plasmid pSL-nspl4, which 

184 includes the TGEV genomic sequence from nt 15,063-21,504, leading to intermediate 

185 plasmids pSL-nspl4-ExoI, pSL-nspl4-ExoIII, pSL-nspl4-ZF-El, pSL-nspl4-ZF-C and 

186 pSL-nspl4-N7MTase. All intermediate plasmids were digested with Nhel and PacI 

187 restriction enzymes and inserts were cloned into the same sites of pBAC-REP2 (46), 

188 generating the mutant replicons pBAC-REP2-nspl4-ExoI, pBAC-REP2-nspl4-ExoIII, 

189 pB AC-REP2-nsp 14-ZF-F1, pBAC-REP2-nspl4-ZF-C, and pBAC-REP2-nspl4- 

190 N7MTase. Alternatively, inserts digested with Nhel and PacI restriction enzymes were 

191 cloned into the same sites of plasmid PBAC-TGEV-S7.1 (C.M. Sanchez, M. Becares, S. 

192 Zuniga, and L. Enjuanes, unpublished results) leading to pBAC-TGEV-S 7 .i-nspl4-ExoI, 

193 pBAC-TGEV-S 7 .i-nspl4-ExoIII, pBAC-TGEV-S 7 i-nspl4-ZF-H, pBAC-TGEV-S 7 .i- 

194 nspl4-ZF-C, and pBAC-TGEV-S 7 .i-nspl4-N7MTase. All cloning steps were checked 

195 by sequencing of the PCR fragments and cloning junctions. For each mutant sequence, 

196 two independent cDNAs were constructed. 

197 

198 Production of a polyclonal antibody specific for TGEV nspl4. An nspl4 protein 

199 with a 6-His tag fused at its N-terminus was expressed in the baculovirus-insect cell 

200 system, using a baculovirus obtained from the Fei Deng laboratory (Wuhan Institute of 

201 Virology, Wuhan, China). Recombinant nspl4 was purified to near homogeneity by 

202 metal chelate affinity chromatography using Ni-NTA agarose (Sigma-Aldrich, Madrid, 

203 Spain) following standard procedures set up in our laboratory (48). Specific polyclonal 

204 antisera were generated by Biogenes GmbH (Germany) after immunization of rabbits 

205 with purified recombinant nsp 14. 

206 
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207 Transfection and recovery of infectious rTGEVs from cDNA clones. BHK-pAPN-N 

208 or BHK-N cells grown to 90% confluence in 35 mm plates (or 12-well plates for RNA 

209 synthesis experiments) were transfected using 3.3 pg of the corresponding pBAC and 

210 10 pi of Lipofectamine 2000 (Invitrogen) per million cells according to the 

211 manufacturer’s specifications. For recovery of infectious recombinant TGEVs 

212 (rTGEVs) from cDNA infectious clones, BFIK-pAPN-N transfected cells were 

213 trypsinized at 6 h post-trasfection (hpt) and plated over confluent ST monolayers grown 

214 in 35 mm plates. After a 2-day incubation period, the cell supernatants were harvested 

215 (passage 0) (49). 

216 Viral RNA presence at passages 0 and 1 was analyzed by RT-PCR. Total cellular RNA 

217 was extracted using the RNeasy Mini kit (Qiagen) following the manufacturer’s 

218 instructions. DNA was removed by treatment of 7 pg of each purified RNA with 20 U 

219 of DNase I (Roche) for 30 min at 37 °C, and DNA-free RNAs were re-purified using 

220 the RNeasy Mini kit (Qiagen). Reverse transcription was performed with the High 

221 Capacity RNA-to-cDNA™ Kit (Life Technologies) according to the manufacturer’s 

222 instructions. PCRs were performed to analyze genomic RNA (gRNA) by amplifying 

223 nspl4 sequences (using the forward primer nspl3fin_VS and the reverse primer 

224 nspl5prin_RS), and mRNA of the N protein using the forward primer SP (5- 

225 GTGAGTGTAGCGTGGCTATATCTCTTC-3') and the reverse primer N-479RS (5'- 

226 TAGATTGAGAGCGTGACCTTG-3'). 

227 Each rTGEV was cloned by three plaque purification steps and the resultant virus was 

228 grown and titrated as previously described (50). Lysis plaque diameter was determined 

229 by measuring the area of the lysis plaques using the software ImageJ (51). 

230 Full genome sequencing of each rTGEV clone was performed using overlapping PCRs 

231 covering the entire TGEV genome (primer sequences available on request). Both 
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232 strands of each PCR product were sequenced, and contigs were assembled using 

233 SeqMan software from DNASTAR Lasergene® package (SeqMan®. Version 8.0.2 

234 DNASTAR. Madison, WI). 

235 

236 Analysis of viral RNA synthesis and cellular gene expression by quantitative RT- 

237 PCR (RT-qPCR). Total intracellular RNA was extracted at 24 hpt from transfected 

238 BHK-N cells, or at the indicated times post-infection from ST cells infected with 

239 rTGEVs. Total RNA was purified with the RNeasy Mini kit (Qiagen) according to the 

240 manufacturer’s specifications. For the analysis of viral RNA synthesis in transfected 

241 cells, an additional step was used in order to remove DNA from samples. For this 

242 purpose, 7 pg of each RNA was treated with 20 U of DNase I (Roche) for 30 min at 37 

243 °C. DNA-free RNA was re-purified using the RNeasy Mini kit (Qiagen). In all cases, 

244 100 ng of total RNA was used as the template for synthesis of cDNA with random 

245 hexamers using the Fiigh-capacity cDNA transcription kit (Life Technologies) 

246 following the manufacturer’s instructions. 

247 Viral RNA levels were evaluated by RT-qPCR using custom TaqMan assays (Life 

248 Technologies) specific for the TGEV genomic RNA (gRNA) and subgenomic mRNA 7 

249 (mRNA-7), both in positive (+) and negative (-) polarity (Table 2), following standard 

250 procedures set up in our laboratory (52). Cellular gene expression was analyzed using 

251 TaqMan gene expression assays (Applied Biosystems) specific for porcine genes 

252 encoding TNF (Ss03391318_gl), IFN-|3 (Ss03378485_ul), IRF-1 (Ss03388785_ml), 

253 MDA5 (Ss03386373_ul), RIG-I (Ss03381552_ul), or TFG-|3 (Ss03382325_ul), and 

254 porcine 2',5' oligoadenylate synthetase 1 (OAS1) was analyzed by using a custom 

255 TaqMan gene expression assay (53). The |3-glucuronidase (GUSB) gene (TaqMan code 

256 Ss03387751_ul) was used as a reference housekeeping gene, since its expression 
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257 remains constant in both infected and non-infected cells (54). Data were acquired with a 

258 7500 real-time PCR system (Applied Biosystems) and analyzed with 7500 software 

259 v2.0.6. The relative quantifications were performed using the 2~ AAO method (55). All 

260 experiments and data analyses were MIQE compliant (56). 

261 

262 Protein analysis by Western blot. Transfected BHK-N cells and infected ST cells 

263 were collected at the indicated times, and cell lysates were obtained by solubilizing cells 

264 in lx SDS gel-loading sample buffer (57). Cell lysates were resolved by denaturing 

265 electrophoresis in NuPAGE 4-12% Bis-Tris gels with 3-morpholinopropane-l-sulfonic 

266 acid (MOPS) SDS running buffer (Invitrogen) and transferred to polyvinylidene 

267 difluoride (PVDF) membranes ( Immobilon®-P PVDF, Merk-Millipore) employing a 

268 Trans-Blot 8 Turbo™ Transfer System (Bio-Rad), using the manufacturer’s 

269 recommended conditions. Membranes were blocked for 1 h at room temperature (RT) 

270 with 5% skim milk in TBS (20 mM Tris-HCl pH 7.5, 150 mM NaCl) and then 

271 incubated with monoclonal antibodies (mAbs) specific for TGEV N protein (mAb- 

272 3DC10) (58), or (3-Actin (Abeam, 1:10,000). Polyclonal antibodies (pAb) against active 

273 caspase 3 protein (Abeam, 1:1000), or TGEV nsp3 (59) were also used. The blots were 

274 then incubated with a horseradish peroxidase-conjugated secondary antibody diluted in 

275 TBS supplemented with 0.1% Tween 20 (TTBS) and 3% skim milk for 1 h at RT. After 

276 extensive washing with TTBS, the immune complexes were detected using Clarity™ 

277 Western ECL Blotting Substrate (Bio-Rad) and the ChemiDoc XRS + System (Bio-Rad), 

278 according to the manufacturer’s instructions. 

279 Protein amounts were estimated by densitometric analysis using ImageLab 4.1 software 

280 (BioRad). At least three different experiments and appropriate gel exposures were used 
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281 in all cases with similar results. In addition, different exposures of the same experiment 

282 were analyzed to assure that data were obtained from films within the linear range. 

283 

284 Immunofluorescence analysis. ST cells were grown on 12 mm glass coverslips in 

285 DMEM 10% FCS to a confluence of 30-50%, and then cell cultures were synchronized 

286 by serum deprivation (60). Briefly, cells were rinsed with phosphate-buffered saline 

287 (PBS) and incubated with DMEM 0.1% FCS for 48 h, as it had been previously 

288 determined that this incubation time allows for 85-90% of ST cells in GO phase. After 

289 serum starvation, cells were released into cell cycle by incubation in DMEM-10% FCS 

290 for 4 h and subsequently mock-infected or infected at a multiplicity of infection (moi) 

291 of 2 with each rTGEV. At 8 or 16 hpi, cells were washed with PBS, fixed with 4% 

292 paraformaldehyde, permeabilized with cold methanol for 10 min and blocked in PBS 

293 with 10% FCS for 45 min at RT. mAh specific for dsRNA (SCICONS English & 

294 Scientific Consulting Kft, 1:200) and pAb specific for TGEV nspl4 (1:1000, see above) 

295 were used. Bound primary antibodies were detected with Alexa Fluor 488- or 594- 

296 conjugated antibodies specific for mouse or rabbit IgG, respectively (1:500, Invitrogen). 

297 Cell nuclei were stained with 4', 6-diamidino-2-phenylindole (DAPI) (1:200, Sigma). 

298 Confocal microscopy was performed using a Leica SP5 laser scanning microscope, and 

299 images were collected and processed with LAS AF software (Leica, Wetzlar, 

300 Germany). Quantification of the intensity of the fluorescence was done by measuring 

301 the mean grey fluorescence of individual cells from grey-scale projections of 

302 microscopy images. The variance of the intensity was also determined and the 

303 normalized optical density variance was calculated as an indicator of signal dispersion, 

304 as previously described (61). Data represents the average of 30 individual cell 

305 measurements. 
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306 

307 Polyinosinicipolycytidylic acid [poly(I:C)] treatment. ST cells were grown to 

308 confluence in 24-well plates and infected with rTGEVs at an moi of 1. At 12 hpi, 

309 infected cells were transfected with 0.5 pg of poly(I:C) (Sigma-Aldrich) per well using 

310 Lipofectamine 2000 (Invitrogen), according to the manufacturer’s instructions. At 4 hpt, 

311 total intracellular RNA was extracted and analyzed by RT-qPCR, as described above. 

312 Total cell RNA integrity was evaluated with a Bioanalyzer 2100 (Agilent Technologies) 

313 following the manufacturer’s recommendations, and analyzed with 2100 Expert 

314 software (Agilent Technologies). 

315 

316 Statistical analysis. Two-tailed, unpaired Student’s /-tests were used to analyze the 

317 difference in mean values between groups. All results were expressed as mean ± 

318 standard deviation; P values <0.05 were considered significant. 

319 

320 RESULTS 

321 Analysis of RNA synthesis by nspl4-ExoN mutants. CoV nspl4 has been proposed 

322 to play a role in RNA synthesis, although the specific mechanism has not been 

323 determined (17, 20). In addition, there is limited information on the role of 

324 Alphacoronavirus ExoN in the infection context. Therefore, a set of four TGEV ExoN 

325 domain mutants was engineered by reverse genetics using a TGEV-derived replicon 

326 (46). Nspl4 is highly conserved among different CoV genera, and all the previously 

327 described motifs were identified in TGEV nspl4 (Fig. 1A). Two of the mutants were 

328 designed to affect the catalytic residues DE/D within the conserved exonuclease motifs 

329 Exol and Exo III, respectively (Fig. 1 A). These mutations are well characterized in vitro 

330 as abolishing ExoN activity in the case of FlCoV-229E and SARS-CoV (17, 18). The 
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331 recent publication of the SARS-CoV nspl4 protein crystal structure has revealed the 

332 presence of two ZFs in the ExoN domain of the protein, which seem essential for the 

333 function of nspl4 (Fig. 1A) (32). Two mutants affecting ZF1 were engineered, designed 

334 to alter the ZF type but theoretically still allowing zinc coordination. The ZF2 motif was 

335 not modified as it overlaps with the catalytic core, and its modification would most 

336 likely influence nspl4 catalytic activity. An additional TGEV mutant replicon 

337 modifying the N7-MTase domain was designed as a control, as it was previously shown 

338 that this mutation has no effect on ExoN activity (32). The mutant replicons were 

339 transfected into BF1K-N cells (52), and replication and transcription levels were 

340 analyzed at 24 hpt by RT-qPCRs measuring the accumulation of gRNA and mRNA-7, 

341 respectively. RNA levels were compared in all cases with those obtained from a non- 

342 replicative replicon, which served to establish background reference levels of 

343 replication and transcription due to the cytomegalovirus promoter. As expected, the 

344 engineered mutation in the N7-MTase domain did not have any impact on RNA 

345 synthesis (Fig. IB). The Exol and ExoIII mutations led to a modest reduction of 

346 approximately 3 fold that was barely statistically significant in both replication and 

347 transcription when compared to the WT replicon, indicating that these mutants were 

348 competent in RNA synthesis (Fig. IB). Interestingly, the engineered mutations in ZF1 

349 led to two different outcomes. Mutant ZF-H showed a modest reduction in replication in 

350 comparison to WT levels, while the level of mRNA-7 was comparable to that of the 

351 non-replicative replicon, indicating that the ZF-H mutation completely abolished 

352 transcription (Fig. IB). In contrast, the ZF-C mutant was competent in both replication 

353 and transcription (Fig. IB). 

354 
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355 Mutations abolishing ExoN activity resulted in a lethal phenotype. To study the 

356 role of TGEV ExoN activity during infection, the Exol and ExoIII mutant infectious 

357 viral cDNAs were generated. The mutant cDNAs were transfected into BHK-pAPN 

358 cells expressing TGEV N protein that were plated over ST cells, susceptible to viral 

359 infection. At 48 hpt, cell supernatants were harvested (passage 0) and used to infect ST 

360 cells (passage 1), following the standard protocol set up in the laboratory (49). No 

361 cytopathic effect was observed after Exol or ExoIII mutant cDNA transfection (data not 

362 shown). In order to analyze whether viral RNA synthesis occurred, total RNA was 

363 extracted from passage 0 and passage 1 and the presence of viral RNA was evaluated by 

364 RT-PCR. Viral genomic RNA (gRNA) corresponding to nspl4 sequences, as well as 

365 the mRNA of N protein (mRNA-N) produced by WT virus were detected both at 

366 passage 0 and 1, indicating an efficient virus recovery (Fig. 2A). In contrast, gRNA and 

367 mRNA-N synthesized from rTGEV ExoN mutants were only amplified from passage 0 

368 samples and not detected at passage 1, indicating a failure in mutant virus recovery (Fig. 

369 2A). Sequencing of the nsp!4 gRNA fragment confirmed that each rTGEV carried the 

370 desired engineered mutations. 

371 In order to determine whether a defect in RNA synthesis prevented recovery of mutant 

372 viruses, RNA accumulation was analyzed by RT-qPCR 24 hours after transfection of 

373 the ExoN mutant cDNAs into BF1K-N cells. Results indicated a decrease in replication 

374 (15 fold reduction) and transcription (2 fold reduction, not statistically significant), 

375 compared to WT cDNAs (Fig. 2B). A similar reduction in the accumulation of viral 

376 RNA was observed both for the positive and negative strands (Fig. 2B). These data are 

377 in agreement with the decrease in viral RNA synthesis observed for FlCoV-229E, 

378 MHV, and SARS-CoV (17,21). Nevertheless, in the case of MHV and SARS-CoV, this 

379 reduction in viral RNA synthesis did not prevent virus recovery. Moreover, previous 
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380 rTGEV mutants obtained in our laboratory that had a higher reduction in RNA synthesis 

381 led to the recovery of viable viruses (62), suggesting that there is another cause for the 

382 lethal phenotype of TGEV Exol and ExoIII mutants. 

383 As nspl4 has also a role in cap synthesis, mutations could have an effect on viral 

384 protein accumulation, although it has been demonstrated that SARS-CoV Exol and 

385 ExoIII mutants were competent in N7-MTase activity (19). Variation in protein 

386 synthesis of the rTGEV mutants was evaluated by Western blot. Accumulation of nsp3, 

387 which is directly translated from gRNA, did not vary significantly with respect to WT 

388 virus despite the reduction of gRNA levels (Fig. 2C), suggesting that viral protein 

389 synthesis was not affected in the Exol and ExoIII mutants. Unfortunately, detection of 

390 other viral structural proteins after cDNA transfection failed even for WT virus, due to 

391 low BAG transfection efficiency. 

392 Altogether, the results indicate that ExoN activity is required for TGEV viability. 

393 Similar results were obtained with HCoV229E (17), therefore the requirement of ExoN 

394 activity may be common to all members of the Alphacoronavirus genus. 

395 

396 Effect of mutations in nspl4 zinc finger 1 in viral phenotype. In order to study the 

397 functional role of nspl4 ZF1 in the context of CoV infection, rTGEV infectious clones 

398 harboring the ZF-H and ZF-C mutations were engineered. As expected, the ZF-H virus 

399 was not recovered (data not shown), probably due to similar viral RNA transcription 

400 defects as those observed with the replicon system (Fig. IB). Interestingly, MHV 

401 substituting two Zn coordinating residues by alanine also led to a lethal phenotype (63), 

402 suggesting an important role for ZF1 in nspl4 function. 

403 The rTGEV-ZF-C virus was successfully recovered with a titer of 4.3x10 7 pfu/ml, 

404 similar to that of the WT virus (8.7 xlO 7 pfu/ml). The introduced mutations were present 
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405 in the recovered mutant virus, even after 10 passages in cell cultures of plaque-purified 

406 viruses, indicating that they were stably maintained in the rTGEV genome. Moreover, 

407 the full-length genomes of two independent clones of rTGEV-WT and rTGEV-ZF-C 

408 were sequenced, and the only differences found between the rTGEV-WT and the 

409 rTGEV-ZF-C were the desired engineered mutations. 

410 Infection with rTGEV-ZF-C induced a delayed cytopathic effect (CPE) in comparison 

411 to rTGEV-WT, and lysis plaques formed by rTGEV-ZF-C were significantly smaller 

412 (half size of diameter) than those from rTGEV-WT (Fig. 3A). This plaque size 

413 reduction was observed both when analyzing uncloned viruses obtained from 

414 transfection and in plaques formed by six individual clones (data not shown). Viral 

415 growth was analyzed at different times post-infection, at both low (0.05) and high (5) 

416 moi. Infectious titer of rTGEV-ZF-C was up to 10 fold higher than that of rTGEV-WT 

417 at early times post-infection, although maximum titer observed from 24 hpi was 

418 identical to that of the rTGEV-WT (Fig. 3B). 

419 Both rTGEV-ZF-C and rTGEV-WT showed the same gRNA accumulation kinetics at 

420 early times post-infection, corresponding with the active RNA synthesis stage of the 

421 TGEV infectious cycle (Fig. 3C), in agreement with our observations using TGEV 

422 replicons (Fig. IB). Nevertheless, at late times post-infection, rTGEV-ZF-C 

423 accumulated lower levels of gRNA than rTGEV-WT (Fig. 3C). This data, together with 

424 the reduced CPE and smaller plaque phenotype, could be attributed to a less efficient 

425 dissemination of the rTGEV-ZF-C in comparison to rTGEV-WT. 

426 

427 Antiviral response elicited by rTGEV-ZF-C. Previous results from our laboratory 

428 indicated that CPE produced by TGEV and lysis plaque size could correlate with 

429 apoptosis induction as a consequence of the host antiviral response (53). TGEV induces 
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430 apoptosis following a caspase-dependent pathway that involves the processing of two 

431 initiator proteases (caspase 8 and 9), as well as three downstream effector caspases 

432 (caspases 3, 6 and 7) (64, 65). To evaluate whether rTGEV-ZF-C CPE and lysis plaque 

433 size reduction were due to reduced apoptosis, caspase 3 activation was evaluated by 

434 Western blot using specific antibodies. At 24 hpi, when rTGEV-ZF-C and rTGEV-WT 

435 infectious titers were similar (Fig. 3B), rTGEV-ZF-C was found to induce lower levels 

436 of active caspase 3 in comparison to the rTGEV-WT (Fig. 4A and B, left pannels). It 

437 has been described that TGEV N protein is cleaved by caspases within the host cell 

438 (66). In agreement with the reduced caspase 3 activation, decreased cleavage of N 

439 protein by caspases was observed during rTGEV-ZF-C infection in comparison to 

440 rTGEV-WT (Fig. 4A and B, right pannels). Altogether, these results indicated that the 

441 rTGEV-ZF-C virus triggered a reduced apoptosis. 

442 TGEV-induced apoptosis is linked to the antiviral response produced during infection 

443 (53, 54). Therefore, accumulation of IFN-(3 and TNF mRNAs was measured as a marker 

444 of the activation of the IRF-3 and NFk| 3 pathways, respectively, during the dsRNA- 

445 triggered antiviral response. rTGEV-ZF-C infection led to significantly reduced levels 

446 of IFN-J3 (up to 25 fold) and TNF (up to 17 fold) in comparison to rTGEV-WT infection 

447 (Fig. 5A). The induction of interferon-stimulated genes (ISGs) was also evaluated, by 

448 measuring the accumulation of IRF-1, OAS, RIG-I, and MDA-5 mRNAs. In agreement 

449 with the reduced IFN-B production, a significant decrease in the expression of these 

450 genes was observed in rTGEV-ZF-C when compared to rTGEV-WT-infected cells (Fig. 

451 5B). In contrast, the accumulation of TGF(5 mRNA, whose expression did not change 

452 with TGEV infection (54), was similar in both cases (Fig. 5B), discarding a general 

453 decrease of host cell transcription after rTGEV-ZF-C infection. 
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454 Overall, these results indicated that rTGEV-ZF-C triggered a reduced antiviral response 

455 in comparison to that induced by rTGEV-WT. 

456 

457 dsRNA accumulation in rTGEV-ZF-C infected cells. In non-immune cells, RIG-I 

458 and MDA-5 cytoplasmic sensors mainly mediate the production of cytokines and IFN. 

459 These proteins recognize the presence of dsRNAs in the cytoplasm, activating 

460 transcription factors IRF-3 and NF-kB, which finally leads to the expression of type I 

461 IFN and proinflammatory cytokines, respectively (67). Because rTGEV-ZF-C infection 

462 led to reduced IFN-J3 and TNF accumulation, dsRNA levels were analyzed by 

463 immunofluorescence in synchronized and infected ST cells. Nspl4 was also detected by 

464 immunofluorescence to analyze the effect of ZF1 mutation on nspl4 protein subcellular 

465 location. The nspl4-ZF-C mutant protein localized similarly to the WT protein (Fig. 

466 6A), with a pattern similar to that of other viral proteins present in the CoV RTC such 

467 as nsp3 (data not shown) (59). 

468 At early times post-infection, no differences in dsRNA distribution pattern were 

469 observed in rTGEV-ZF-C and rTGEV-WT infected cells (Fig. 6A, 8 hpi panels). Co- 

470 localization between nspl4 and dsRNA was observed, probably corresponding to active 

471 sites of viral RNA synthesis (Fig. 6A, lower 8 hpi panels). At late times post-infection, 

472 dsRNA labeling produced a scattered pattern in rTGEV-WT infected cells, uniformly 

473 distributed throughout the cytoplasm (Fig. 6A, middle 16 hpi panels). This dsRNA 

474 pattern was previously observed in CoV infection at late hpi, (10, 68). In contrast, 

475 dsRNA was still concentrated in well-defined perinuclear granules in rTGEV-ZF-C 

476 infected cells (Fig. 6A, 16 hpi panels). In order to quantify the cytoplasmic dsRNA 

477 accumulation and distribution, the mean intensity of fluorescence as well as the 

478 normalized variance of the intensity of fluorescence were measured in thirty individual 
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479 cells. In the case of a uniform distribution of the signal the variance is minimal, while 

480 its value increases when the signal is irregularly distributed (61). An un-biased analysis 

481 showed that, at early times post-infection, no significant differences were observed in 

482 both mean and variance of intensity between rTGEV-ZF-C and rTGEV-WT infected 

483 cells (Fig. 6B). In contrast, significant reduction in the intensity of fluorescence was 

484 observed in rTGEV-ZF-C infected cells compared to the rTGEV-WT (Fig 6B). 

485 Moreover, the quantification of the signal distribution resulted in a significant increase 

486 in the normalized variance in the case of rTGEV-ZF-C infected cells, which meant a 

487 reduction in the spreading of the signal in comparison to that of the rTGEV-WT 

488 infected cells (Fig 6B). This result suggested that the nspl4-ZF-C mutation interfered 

489 with dsRNA accumulation at late times post-infection, which may be the cause for the 

490 reduced induction of subsequent antiviral responses. 

491 

492 Modulation of innate immune response by nspl4-ZF-C. To further investigate 

493 whether rTGEV-ZF-C actively inhibits the cytoplasmic accumulation of dsRNA and, as 

494 a consequence, the activation of the subsequent antiviral response, the dsRNA-sensing 

495 pathway was activated by the addition of exogenous dsRNA in the context of viral 

496 infection. 

497 In order to evaluate the potential of rTGEV-WT or rTGEV-ZF-C viruses in the 

498 inhibition of dsRNA-triggered antiviral responses, ST cells were infected with rTGEV- 

499 ZF-C or rTGEV-WT at an moi of 1 and at 12 hpi transfected with the dsRNA analog 

500 poly(I:C). At 4 hpt, viral RNA accumulation and the induction of the innate immune 

501 response were analyzed by measuring gRNA, IFN-J3, and TNF mRNA accumulation by 

502 RT-qPCR. In the absence of poly(I:C), IFN-J3 and TNF mRNA production were 

503 significantly reduced in rTGEV-ZF-C infected cells, as described above (Fig. 7A). 
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504 Poly(I:C) treatment did not significantly affect virus production, as similar gRNA levels 

505 were detected in the absence and presence of poly(I:C) (Fig. 7A). rTGEV-WT virus 

506 together with poly(I:C) triggered an exacerbated innate immune response, as expected 

507 from the combined action of rTGEV-WT infection and the dsRNA analog poly(I:C) 

508 (Fig. 7A). In contrast, this response was significantly reduced in the case of rTGEV-ZF- 

509 C infection, despite the exogenous activation by poly(I:C) (Fig. 7A). 

510 To test whether during rTGEV-ZF-C mutant virus infection the decreased dsRNA levels 

511 led to a reduced recognition and, subsequently to a diminished innate immune response, 

512 RNase L activation was analyzed. As previously described, rTGEV-WT infection 

513 hardly produces RNase L activation and RNA degradation (53). In agreement with this 

514 data, no significant difference in RNA degradation was observed for the rTGEV-ZF-C 

515 mutant infected cells compared with those infected with the parental virus (Fig. 7B). In 

516 contrast, when RNase L activation was triggered by exogenous dsRNA [poly(I:C)] 

517 transfection together with infection, a significant decrease was observed in RNA 

518 degradation during rTGEV-ZF-C infection compared with rTGEV-WT infection (Fig. 

519 7B). These results indicated that indeed during rTGEV-ZF-C infection there was a 

520 decreased dsRNA recognition, as expected due to the lower dsRNA accumulation, 

521 which, in consequence, caused a reduced IFN production. 

522 In a complementary approach, co-infections with rTGEV-ZF-C and rTGEV-WT were 

523 performed in order to evaluate whether the antiviral response elicited by rTGEV-WT 

524 infection was modulated by the rTGEV-ZF-C virus. ST cells were infected at different 

525 moi of each virus, and IFN-J3 and TNF mRNA accumulation were evaluated by RT- 

526 qPCR at 16 hpi. A slight difference in gRNA levels between rTGEV-ZF-C and rTGEV- 

527 WT, which as shown above does not influence infectious titers, was also observed 

528 during independent infections (Fig. 8). No significant differences in gRNA 
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529 accumulation were detected during co-infections (Fig. 8). Moreover, the nspl4 region 

530 was sequenced and the proportion of each rTGEV at the moment of the analysis was 

531 estimated to be similar to the one in the inoculum (data not shown). A decreased innate 

532 immune response was observed in cells co-infected with rTGEV-ZF-C and rTGEV-WT 

533 at a 1:1 ratio, and more significantly in co-infections at a 5:1 ratio, when compared to 

534 the same dose of rTGEV-WT alone (Fig. 8). This result suggests that the nspl4-ZF-C 

535 mutation could be acting in trans, reducing the dsRNA produced during rTGEV-WT 

536 infection and, as a consequence, decreasing the subsequent innate immune response. 

537 Furthermore, this data strongly suggests that the mutant virus actively reduced the 

538 antiviral response triggered by the parental virus, as IFN13 and TNF accumulation were 

539 not the sum of those induced by each virus independently and, in fact, the reduction was 

540 dependent on the moi used for the rTGEV-ZF-C virus. 

541 Altogether these data indicate that rTGEV-ZF-C leads to a reduction in accumulation of 

542 dsRNA, either exogenous or that produced during infection. Our results indicate that 

543 CoV nspl4 modulates the innate immune response most likely by affecting dsRNA 

544 accumulation. 

545 

546 DISCUSSION 

547 Previous studies have shown that CoV ExoN activity is required for the high-fidelity 

548 replication of CoV genome (21, 25). Although there is extensive knowledge on 

549 Betacoronavirus ExoN functions during infection, the function of Alphacoronavirus 

550 ExoN during infection has not been analyzed apart from biochemical characterization in 

551 vitro. In this study, we generated a set of mutants in different domains of TGEV ExoN, 

552 based on the published structure of npsl4 (32). Mutants with abrogated ExoN activity 

553 were lethal, despite being competent in viral RNA and protein synthesis. The reduction 
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554 in viral RNA synthesis reported here for TGEV and previously for HCoV229E (17) was 

555 similar to that observed for the analogous equivalent viable mutants in Betacoronavirus, 

556 such as SARS-CoV (21). Therefore, the reduction in RNA synthesis is most likely not 

557 itself the causative factor for the lethality of the ExoN mutations. 

558 Mutagenesis of a residue in the ZF1 motif, included in the ExoN domain of nspl4, led 

559 to viable viruses reaching an infectious titer similar to that of the parental virus. The 

560 rTGEV-ZF-C mutant accumulated lower amounts of dsRNA in the cytoplasm, and 

561 subsequently triggered a reduced antiviral response, as evidenced by decreased levels of 

562 IFN-J3 and TNF mRNAs. This work showed for the first time that CoV ExoN domain, 

563 beyond its well-known role in proofreading, participates in innate immune response 

564 modulation. 

565 The relevance of nspl4 ZF1 motif was previously suggested by biochemical studies 

566 (32), as well as reverse genetics analyses, which led to non-recovery of a virus carrying 

567 mutations to alanine of the ZF1 domain residues responsible for zinc coordination (63). 

568 Structural data, together with the fact that those mutations severely altering the ZF were 

569 lethal in MF1V, whereas MF1V mutants lacking ExoN activity were viable (21, 63), 

570 strongly suggest that this domain is not directly involved in ExoN catalytic activity. 

571 Although the ZF1 motif may have a role in the maintenance of specific interactions of 

572 the nspl4 protein with its RNA substrate or other RTC proteins. Our data support the 

573 relevance of this domain, as conservative mutations theoretically allowing zinc 

574 coordination, such as the ZF-F1 mutant, completely abolished transcription in both 

575 replicon and virus systems. In addition, our data were also in agreement with ZF1 motif 

576 not being directly involved in ExoN catalytic activity, as preliminary sequencing data 

577 has shown similar mutation rates between rTGEV-WT and rTGEV-ZF-C viruses (data 

578 not shown). 
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579 As other ZF motifs, the ExoN ZF1 domain may have a role in RNA-nspl4 interaction, 

580 and therefore mutations in ZF1 may impact its properties in terms of affinity and/or 

581 specificity. In vitro studies have revealed that CoV nspl4 binds to both ssRNA and 

582 dsRNA with no sequence specificity, although data indicated that ExoN activity is much 

583 higher in substrates with dsRNA structure (22). ZFs are typically involved in protein 

584 binding to nucleic acids (35), though unfortunately the recent nspl4 crystallographic 

585 structure does not include an RNA substrate, so the specific RNA-protein interactions 

586 involving each ZF motif were not elucidated (32). The molecular bases for the 

587 mechanism by which ZF proteins recognize dsRNA are complex, as they are highly 

588 dependent on contact with phosphates and hydrophobic interactions (69). This 

589 complexity makes it difficult to predict how mutations may impact the binding 

590 properties of a ZF-containing protein to its substrate. Further analyses on the nspl4 ZF- 

591 C mutant protein exonuclease or N7-MTase activities, or RNA binding ability, are 

592 required to decipher the role of nspl4 in the dsRNA response modulation. This would 

593 require in vitro analysis with the purified WT and ZF-C nspl4 proteins (17, 28). 

594 Moreover, to determine the role on ExoN activity purified nsplO should also be 

595 included (22). This would represent an extensive biochemical characterization that 

596 would be the aim of future work. 

597 In spite of reaching similar infectious titer and gRNA accumulation, cells infected with 

598 mutant rTGEV-ZF-C accumulated lower levels of dsRNA than the rTGEV-WT infected 

599 cells. Two alternative explanations are compatible with this fact: (i) the ZF-C mutant 

600 degrades dsRNA more efficiently than the WT virus, or (ii) the ZF-C mutant generates 

601 lower amounts of dsRNA intermediates during infection. Nevertheless, the biochemical 

602 activity of nspl4 as an RNase, which has been clearly reported in vitro (17, 22), is 

603 hardly compatible with a role in generating dsRNA. Moreover, the fact that rTGEV-ZF- 
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604 C reduced the antiviral response triggered both by exogenous dsRNA [poly(I:C)] or 

605 infection-produced dsRNA (during co-infections of WT and mutant viruses) makes the 

606 first hypothesis more feasible. Therefore, our data suggests that nspl4 ExoN plays a 

607 role in degrading dsRNA intermediates that could act as PAMPs, triggering the 

608 subsequent immune response. In this working model (Fig. 9), the engineered rTGEV- 

609 ZF-C mutant was more efficient in this task, leading to reduced amounts of dsRNA and, 

610 subsequently, inducing a very weak innate immune response. 

611 Theoretically, it would be expected that a reduced antiviral response should lead to an 

612 increase in viral replication, eventually leading to higher virus titers. On the other hand, 

613 as viruses evolve leading to optimum replication in a cell host, it is quite frequent that, 

614 in cell cultures any disturbance in host antiviral response does not cause a significant 

615 effect in virus titers. This seems to be the case for TGEV virus infecting porcine ST 

616 cells, as mutant viruses causing increased antiviral response replicate to similar titers as 

617 wild-type virus in tissue cultures (53). In line with that observation, the decreased IFN 

618 production caused by the rTGEV-ZF-C mutant virus did not lead to an increased viral 

619 titer. It is worth noting that rTGEV-WT infection of ST cells produce around 35 pg/ml 

620 of IFN13 (54). This concentration is much lower than that required to affect TGEV 

621 replication. In fact, we have found that a concentration of 9000 pg/ml of IFNB is 

622 required to cause a 20-fold decrease in viral titer (M. Becares, S. Zuniga and L. 

623 Enjuanes, unpublished data). Therefore, it was not surprising that in our experimental 

624 conditions, it was not possible to detect an effect on virus titers of the ZF-C mutation. 

625 To our knowledge, this work represents the first link of CoV ExoN activity and the 

626 degradation of dsRNA leading to innate immune suppression, although this function of 

627 CoV nspl4 ExoN was recently hypothesized (42). Beyond CoVs, only arenaviruses 

628 encode a protein with exonuclease activity within their genome. Similar to CoV nspl4, 
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629 the arenavirus nucleoprotein (NP) has cap binding function and carries 3'-5' 

630 exonuclease activity, degrading short RNA molecules and inhibiting IFN production in 

631 arenavirus infected cells. Studies revealed that mutants lacking the 3'-5' exonuclease 

632 activity do not inhibit IFN induction, strongly implicating it in modulation of cellular 

633 antiviral responses (70-72). Additionally, it has been shown that the arenavirus NP 

634 directly interacts with the cytosolic RNA sensors RIG-I and MDA-5 (73), giving 

635 support to the hypothesis that degradation of dsRNA PAMPs by the arenavirus 

636 exonuclease avoids further recognition by RIG-I or MDA5. In the case of CoVs, which 

637 efficiently inhibit IFN induction during infection, nspl4 could also be acting by 

638 degrading dsRNA PAMPs. Interestingly, CoV nspl4 is known to interact with the 

639 cellular helicase DDX1, which has recently been described as a cytosolic dsRNA sensor 

640 that activates type I IFN responses (74). 

641 In agreement with our working model hypothesis (Fig. 9), CoV mutants lacking ExoN 

642 activity would trigger an exacerbated innate immune response. An enhanced antiviral 

643 response, together with the reported defects in viral RNA synthesis and the increased 

644 mutation rate, could contribute to the lethal phenotype observed for Alphacoronaviruses 

645 lacking ExoN activity. In fact, SARS-CoV ExoN mutants were rescued in Vero cells, 

646 defective in IFN synthesis, and MHV ExoN mutants rescue was initially impaired when 

647 IFN competent cells were used (63, and M. Denison, personal communication). SARS- 

648 CoV lacking ExoN activity are attenuated in young and aged diseased mice, showing a 

649 faster clearance than the WT virus (75). This phenotype could be compatible with the 

650 triggering of a higher innate immune response by these mutants, although this aspect 

651 was not analyzed. 

652 Further studies are needed to understand the role of ribonucleases, both ExoN and 

653 NendoU, in antagonizing the activation of the dsRNA-induced antiviral response in the 
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654 different CoV genera. These proteins may counteract the innate immune response by 

655 degrading dsRNA and limiting activation of antiviral pathways by cytosolic sensors, 

656 and therefore they may serve as targets for the design of novel antiviral strategies. 

657 
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901 FIGURE LEGENDS 

902 Figure 1. RNA synthesis of TGEV nspl4 mutants. (A) Schematic representation of 

903 TGEV nspl4. Exonuclease (ExoN) andN7 methyltransferase (N7-MTase) domains are 

904 indicated. Motifs I, II and III conferring the ExoN active site are shown (grey boxes). 

905 The three zinc fingers (ZF1, ZF2 and ZF3) are indicated (black boxes), as the amino 

906 acids forming the S-adenosylmethionine (SAM) binding pocket in the N7-MTase 

907 domain (white box). Within the partial sequence alignments, positions of key amino 

908 acids in each motif (white letters in balck boxes) are indicated. Those mutated in each 

909 rTGEV construct with respect to the wild-type (WT) sequence are shown as black 

910 letters in grey boxes. (B) Quantification of replication and transcription levels in each 

911 rTGEV mutant replicon. Non-replicative (NR) and wild-type (WT) replicons were 

912 included as controls. Genomic RNA (gRNA) and subgenomic mRNA from gene 7 

913 (mRNA-7) were analyzed by RT-qPCR using specific TaqMan assays after transfection 

914 of each replicon in BF1K-N cells. Mean values from eight independent transfection 

915 experiments are plotted; error bars represent standard deviation. *, p-value < 0.05 ; **, 

916 p-value < 0.01, ***, p-value<0.001. 

917 

918 Figure 2. Analysis of rTGEV ExoN mutants. (A) RT-PCR analysis of viral RNA 

919 from passage 0 and passage 1 of rTGEV-ExoI and rTGEV-ExoIII mutants, compared 

920 with rTGEV-WT. Genomic (gRNA) and subgenomic mRNA of gene N (mRNA-N) 

921 were analyzed. A and B indicate duplicate infectious cDNAs, tested individually for 

922 each mutant. Size (bp) of molecular weight markers (Mw) is shown on the right. (B) 

923 Quantification of genomic RNA (gRNA) and mRNA of gene 7 (mRNA-7). Both 

924 negative and positive strands were measured by RT-qPCR using specific TaqMan 

925 assays 24 h post-transfection of infectious cDNAs into BF1K-N cells. Means from four 
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926 independent transfection experiments are plotted; error bars represent standard 

927 deviation. (C) Detection of viral proteins by Western blot. Total protein was extracted 

928 24h after transfection of BHK-N cells with infectious cDNAs from Exol, ExoIII, WT, a 

929 non-replicative rTGEV (NR), or mock transfected cells (M). Viral protein nsp3 and (3- 

930 actin (as loading control) were detected using specific antibodies. Right panel represents 

931 the quantification of the bands by densitometry, corrected by the amount of |3-actin. 

932 Means from two independent transfection experiments are plotted; error bars represent 

933 standard deviation. *, p-value < 0.05 ; ***, p-value < 0.001. 

934 

935 Figure 3. Growth of nspl4 mutant rTGEV-ZF-C in tissue cultures. (A) Lysis 

936 plaques produced by mutant rTGEV-ZF-C and rTGEV-WT at 48 h post-infection in ST 

937 cells. Right panel represents the mean diameter of 10 individual lysis plaques; error bars 

938 indicate standard deviation. (B) Growth kinetics of rTGEV-WT and rTGEV-ZF-C 

939 viruses. ST cells were infected at low (0.05, left) and high (5, right) moi with mutant 

940 rTGEV-ZF-C (ZF-C) or rTGEV-WT (WT); supernatants were collected at different 

941 times post-infection, and infectious titers were determined by plaque titration on ST 

942 cells. Means from three independent experiments are plotted; error bars represent 

943 standard deviation. (C) RNA synthesis of rTGEV-WT and rTGEV-ZF-C viruses. ST 

944 cells were infected at an moi of 1 with mutant rTGEV-ZF-C or rTGEV-WT; 

945 intracellular RNA was collected at different times post-infection, and genomic RNA 

946 (gRNA) was quantified by RT-qPCR using specific TaqMan assays. Means from three 

947 independent experiments are plotted; error bars represent standard deviation. *, p-value 

948 <0.05 ; **, p-value<0.01; ***, p-value < 0.001. 

949 
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950 Figure 4. Induction of apoptosis by mutant rTGEV-ZF-C. (A) Detection of active 

951 caspase 3 (casp3) and TGEV N protein (N) by Western blot. Total protein was extracted 

952 from ST cells 24 h after infection with mutant rTGEV-ZF-C, rTGEV-WT, or mock- 

953 infected (M). A and B indicate duplicate viral clones that were tested individually. 

954 Casp3, N, and (3-actin (loading control) were detected using specific antibodies; 

955 procaspase 3 (procasp3) and the form of TGEV N protein cleaved by caspases (N-cl) 

956 are also indicated. (B) Quantification of the bands by densitometry, corrected by 

957 amount of the (3-actin. Relative levels of protein (r.u.) were based on comparison with 

958 the WT virus, which was considered to be 100%. Means from three independent 

959 experiments are plotted; error bars represent standard deviation. *, p-value < 0.05 ; **, 

960 p-value<0.01. 

961 

962 Figure 5. Innate immune response induced by mutant rTGEV-ZF-C. ST cells were 

963 infected with rTGEV-ZF-C or rTGEV-WT virus at an moi of 1, and intracellular RNA 

964 was collected at different times post-infection. (A) Quantification of IFN-(3 and TNF 

965 mRNA was performed by RT-qPCR using specific TaqMan assays; relative mRNA 

966 levels were based on comparison with mock-infected cells. Means from three 

967 independent experiments are plotted; error bars represent standard deviation. (B) 

968 Induction of interferon-stimulated genes by rTGEV-ZF-C. Quantification of IRF-1, 

969 OAS, RIG-I, MDA5 and TGFB mRNAs was performed by RT-qPCR using specific 

970 TaqMan assays; relative mRNA levels were based on comparison with mock-infected 

971 cells. Numbers below asterisks indicate the fold change in induction by rTGEV-WT 

972 relative to that of mutant rTGEV- ZF-C. Means from four independent experiments are 

973 plotted; error bars represent standard deviation. *, p-value < 0.05 ; **, p-value<0.01; 

974 ***, p-value < 0.001. 
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975 

976 Figure 6. Accumulation of dsRNA in cells infected with mutant rTGEV-ZF-C. 

977 Confocal microscopy analysis was performed on synchronized ST cells infected with 

978 the rTGEV-ZF-C or rTGEV-WT virus at an moi of 1. (A) Immunofluorescence images 

979 of ST cells mock-infected (mock) or infected with rTGEV-ZF-C (ZF-C) or rTGEV- WT 

980 at 8 and 16 h post-infection. Nspl4 was detected using a specific polyclonal antisera 

981 and a secondary antibody staining green; dsRNA was detected using monoclonal 

982 antibody mAb-J2 and a secondary antibody staining red; DAPI (4',6'-diamidino-2- 

983 phenylindole) (blue) was used to stain the nuclear DNA. Co-localization is indicated by 

984 yellow pixels in the merge panels. (B) Quantification of the mean intensity of 

985 fluorescence and normalized standard deviation of fluorescence intensity (intensity 

986 variance) of nspl4 and dsRNA. Means of 30 individual cells plotted; error bars 

987 represent standard deviation. ***, p-value < 0.001 

988 

989 Figure 7. Modulation of dsRNA-induced antiviral response by mutant rTGEV-ZF- 

990 C. (A) ST cells were mock-infected (white bars) or infected with mutant rTGEV-ZF-C 

991 (grey bars) or rTGEV-WT (black bars) at an moi of 1, and subsequently transfected at 

992 12 h post-infection with poly(I:C). At 4 h post-transfection, total intracellular RNA was 

993 collected, and quantification of viral genomic RNA (gRNA), IFN-J3, and TNF mRNAs 

994 was performed by RT-qPCR using specific TaqMan assays. Relative mRNA levels 

995 were based on comparison with mock-infected, non-transfected cells. Means from three 

996 independent experiments are plotted; error bars represent standard deviation. (B) 

997 Cellular RNA integrity. Total RNA was extracted from mock-infected or rTGEV-WT 

998 (WT) and rTGEV-ZF-C (ZF) infected ST cells, non-treated (-PolylC) or treated 

999 (+PolyIC) with poly(I:C). The RNA was then analyzed using a Bioanalyzer. The 
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1000 position of 28S and 18S rRNAs is indicated (left panel). Graph of 28S rRNA integrity, 

1001 as measured by the Bioanalyzer. Error bars indicate the standard deviation from four 

1002 independent experiments, f.u., fluorescence units. *, p-value < 0.05; **, p-value<0.01; 

1003 ***, p-value<0.001. 

1004 

1005 Figure 8. Modulation of rTGEV-WT induced antiviral response by rTGEV-ZF-C 

1006 mutant. ST cells were mock-infected, or individually infected with mutant rTGEV-ZF- 

1007 C or rTGEV-WT, or co-infected with each virus at two different ratios (1:1 or 5:1 ZF- 

1008 C:WT ratio; the moi of each virus used is indicated by numbers in X-axis). At 16 h 

1009 post-infection, total intracellular RNA was collected, and quantification of viral 

1010 genomic RNA (gRNA), IFN-J3, and TNF mRNAs was performed by RT-qPCR using 

1011 specific TaqMan assays. Relative mRNA levels were based on comparison with mock- 

1012 infected cells. Means from three independent experiments are plotted; error bars 

1013 represent standard deviation. *, p-value < 0.05; **, p-value<0.01. 

1014 

1015 Figure 9. Working model for the role of nspl4 in the counteraction of antiviral 

1016 responses during CoV infection. Schematic overview of the host cell dsRNA-induced 

1017 antiviral pathway. The proposed mechanism of action for CoV nspl4 is indicated. 

1018 
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Table 1. Primers used for nsp!4 mutagenesis 


MUTANT 

PRIMER a) 

SEQUENCE (5’-3’) (b) 

nspl4-ExoI 

ExoIVS 

ExoIRS 

CTTGGGTTTGCAGTTGCTGGTGCACATG 

CATGTGCACCAGCAACTGCAAACCCAAG 

nspl4-ExoIII 

ExoIIIVS 

ExoIIIRS 

GCTAGTGGTGCGGCTATCATGACTAG 

CTAGTCATGATAGCCGCACCACTAGC 

nspl4-ZF-H 

CHCHVS 

CHCH RS 

CAAAAATGTGAACACGGCAAAAG 

CTTTTGCCGTGTTCACATTTTTG 

nspl4-ZF-C 

CCCC_VS 

CCCCRS 

GCTTGCTTCAAGTGTGCATTAGGATG 

CATCCTAATGCACACTTGAAGCAAGC 

nsp!4-N7MTase 

N7MTase_VS 

N7MTase_RS 

GCTTGCTTCAAGTGTGCATTAGGATG 

GATTACCCACCGCGTGAATTGCAGC 


(a> VS, forward primer; RS, reverse primer. 


(b) Mutant nucleotides are in bold. 
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Table 2. TaqMan assays used for RT-qPCR of viral gRNA and sgmRNA 


AMPLICON 

TAQMAN ASSAY a 

SEQUENCE (5’-3’) <b ' 


PrimerVS 

Rep VS 

TTCTTTTGACAAAACATACGGTGAA 

gRNA(+) 

PrimerRS 

RepRS 

CTAGGCAACTGGTTTGTAACATCTTT 


Probe 

Rep-MGB 

FAM-AGGGCACCGTTGTCA-MGB 


PrimerVS 

RepRS 

CTAGGCAACTGGTTTGTAACATCTTT 

gRNA(-) 

PrimerRS 

Rep VS 

CTAGGCAACTGGTTTGTAACATCTTT 


Probe 

cRep-MGB 

FAM-CTGTTCACCGTATGTTT-MGB 


PrimerVS 

Ldrt-VS 

CGTGGCTATATCTCTTCTTTTACTTTAACTAG 

mRNA7(+) 

PrimerRS 

7 (38)-RS 

AAAACTGTAATAAATACAGCATGGAGGAA 


Probe 

mRNA7-MGB 

FAM-CGAACTAAACGAGATGCT-MGB 


PrimerVS 

7 (38)-RS 

AAAACTGTAATAAATACAGCATGGAGGAA 

mRNA(-) 

PrimerRS 

Ldrt-VS 

CGTGGCTATATCTCTTCTTTTACTTTAACTAG 


Probe 

cmRNA7-MGB 

FAM-AGCATCTCGTTTAGTTCGAGT-MGB 


a VS, forward primer; RS, reverse primer; MGB, minor groove binder group. 
b FAM, 6-carboxyfluorescein. 
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